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We study the temperature dependence of the single particle spectral function as well as of the 
dynamical spin and charge structure factors for the one-dimensional Hubbard model using the finite 
temperature auxiliary field quantum Monte Carlo algorithm. The parameters of our simulations are 
chosen so to at best describe the low temperature photoemission spectra of the organic conductor 
TTF-TCNQ. Defining a magnetic energy scale, Tj, which marks the onset of short ranged 2k / 
magnetic fluctuations, we conclude that for temperatures T < Tj the ground state features of the 
single particle spectral function are apparent in the finite temperature data. Above Tj spectral 
weight transfer over a scale set by the hopping t is observed. In contrast, photoemission data points 
to a lower energy scale below which spectral weight transfer occurs. Discrepancies between Hubbard 
model calculations and experiments are discussed. 

PACS numbers: 71.27. +a, 71.10.-w, 71.10.Fd 



I. INTRODUCTION 



The hallmark of Luttinger liquids lies in spin-charge 
separation; an electron fractionalizes into a spinon car- 
rying the spin degrees of freedom and a holon carrying 
the charge. Detecting spin-charge separation relies on 
the study of the single particle spectral function. From 
the theoretical point of view, the ground state low energy 
properties of the single particle spectral function can be 
obtained from bosonization [l| . This approach yields two 
branch cuts, corresponding to the spinon and the holon, 
dispersing linearly from the Fermi wave vector with spin 
and charge velocities. Beyond this low energy limit, ex- 
act calculations of the spectral function have been carried 
out for the infinite U Hubbard model 2]. Taking into 
account that in this limit the spin velocity vanishes, the 
results stand in agreement with the low energy bosoniza- 
tion picture. Furthermore the calculations reveal higher 
energy features such as a holon shadow band. Hence, 
distinct signatures of Luttinger liquids may be found in 
a wide energy range thus facilitating detection in photoe- 
mission experiments. Beyond the infinite U limit, numer- 
ical simulations such as dynamical DMRG (DDRMG) Q 
or Quantum Monte Carlo 4] can be used to investigate 
the zero temperature properties of the spectral function. 
In particular, T = DDMRG results for the Hubbard 
model have been compared successfully with low temper- 
ature, T = 60K, photoemission experiments on the or- 
ganic one-dimensional conductor TTF-TCNQ |5| . In the 
temperature range 60K < T < 260K experiments point 
towards substantial spectral weight transfer. Keeping the 
model parameters which reproduce the low temperature 
data, our aim is to understand if the experimentally ob- 
served temperature behavior of the spectral function can 
be reproduced by finite temperature model calculations. 



The Hubbard model we consider reads: 

H = -* £ (4,^ + h - c ^ + 

+ U ^n iT n a -/i^(n iT +n u ) (1) 

i i 

where t it the hopping amplitude, U the Coulomb re- 
pulsion, fx the chemical potential, and the first sum runs 
over nearest neighbors. c\ G (c^) creates (annihilates) an 
electron in the Wanier state centered around lattice site 
i and with z-component of spin a =|, j. Comparison 
between DDMRG results and experiments point to a pa- 
rameter set U/t = 4.9, t = OAeV and n = 0.59 for an ad- 
equate description of the TCNQ chain. Throughout this 
article we will keep those parameters fixed and vary the 
temperature. The organization and main results of the 
paper are as follows. In section |H] we briefly present the 
finite temperature auxiliary field quantum Monte Carlo 
(QMC) method and the maximum entropy method we 
have used to analytically continue the imaginary time 
QMC data. Section lllll is dedicated to the results. To 
map out the scales involved in the problem, we first con- 
sider the temperature dependence of the spin and charge 
susceptibilities as well as of the spin and charge dynam- 
ical structure factors. This allows us to define a mag- 
netic crossover energy scale, Tj, below which the mag- 
netic 2k f correlation length increases substantially as a 
function of decreasing temperature. In section UlI Bl we 
analyze the temperature dependence of the single parti- 
cle spectral function and arrive to the conclusion that Tj 
is the only low energy scale at hand in the problem. That 
is, for temperatures below Tj the overall features of the 
zero temperature spectral function are well reproduced. 
Above Tj we observe spectral weight transfer over en- 
ergy scale set by t. With t = OAeV our estimated value 
of Tj is Tj ~ 400K. In the conclusion (Sec. HVJ), we 
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FIG. 1: Comparison of different Maximum Entropy methods 
for the spectral function at k = 0, of the one-dimensional 
Hubbard model at U/t = 4.9, n = 0.59 and (3t = 7. Note the 
semi- log scale. 
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FIG. 2: Spin and charge static susceptibilities, Xs(/c,cj = 0), 
as a function of wave vector. At 2k f the different lines 
from bottom to top correspond to the temperatures (3t = 
2,4,7,10,15,20 



discuss possible scenarios to understand the discrepancy 
between experimental data and model calculations. 



II. QUANTUM MONTE CARLO ALGORITHM 

We have used a generic implementation of the finite 
temperature grand canonical auxiliary field algorithm |(| 
to compute imaginary time displaced Green functions, 
G(fc, r) = (cg(r)ct(O)), as a function of temperature for 
the one-dimensional Hubbard model on a 48-site chain 
with periodic boundary conditions. The spectral function 
is extracted from the imaginary time data by inverting 
the equation: 

(cg(r)ct(O)) = J duoK(r,uj)A(k,uj) with 
1 e~ TU 



K(t,u) 



tt 1 + e-^' 



(2) 



Since this inversion is numerically ill defined we have used 
the Maximum Entropy method and favored a recently 
proposed stochastic version (jj. The stochastic formula- 
tion has the appealing property that it formally contains 
the generic Maximum Entropy method (MEM) |a, |2( at 
the mean field level. In the generic MEM, there is no 
free parameter. In particular, a which determines how 
much information is taken from the default model is de- 
termined self-consistently. In contrast, in the stochastic 
approach, there is no sharp way of determining a and we 
have used the criterion proposed by K. Beach Q. In Fig. 
0we compare the different Maximum Entropy methods 
for the one-dimensional Hubbard model at U/t = 4.9, 
n = 0.59 and /3t = 7. The Brayn and Classic Maximum 
entropy formulations Q yield identical spectral function. 
The two peak structure at uo/t < corresponding to the 
holon and spinon branches, is sharper in the stochastic 
approach. At uo/t > the stochastic spectrum shows less 
features than the Classic Maximum Entropy method. It 
is know that the Classic Maximum Entropy has difficul- 
ties in reproducing flat spectra and generates a curve os- 
cillating smoothly around the correct value. This prob- 
lem seems to be alleviated by the stochastic approach. 
Hence, our overall opinion is that the stochastic approach 
does better at reproducing sharp features as well as flat 
regions in the spectra. Finally we note that we have al- 
ways taken the covariance matrix into account. 



III. RESULTS 

We compute the single particle spectral function as 
well as the dynamical spin and charge structure factors 
as a function of temperature. We choose a parameter set 
which at best describes the low temperature properties 
of the photoemission spectra of TTF-TCNQ |. That is 
for the TCNQ band, U/t = 4.9, t = 0.4 eV and a filling 
fraction n = 0.59. We cover the following range of inverse 
temperatures (3t = 2, 4, 7, 10, 15 and (3t = 20. Below, we 
will first discuss two particle properties so as to pin down 
scales and then consider the temperature behavior of the 
single particle spectral function. 
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FIG. 3: Dynamical spin- spin correlations as a function of tem- 
perature on a logarithmic intensity scale. The inverse tem- 
peratures from top to bottom read: fit = 20, 15, 10, 4. 
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FIG. 4: Dynamical charge-charge correlations as a function 
of temperature on a logarithmic intensity scale. The inverse 
temperatures from top to bottom read: (3t = 20, 15, 10, 4. 
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Dynamical spin and charge correlation 
functions. 



To investigate the spin and charge dynamics we con- 
sider the dynamical susceptibility 
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Fig. [3 plots 

the static spin and charge susceptibilities, Xj(&> = 0) as 
a function of temperature. Those quantities measure cor- 
relation lengths and allow us to identify crossover energy 
scales below which spin and charge fluctuations grow as 
a function of decreasing temperature. Since one dimen- 
sional systems are critical at T = both 2k f spin and 
charge static susceptibilities diverge at T = 0. As appar- 
ent from Fig. 0the crossover scale marking the growth 
of 2k f spin fluctuations is given by Tj ~ tj 10. A sim- 
ilar energy scale for the 2k f charge fluctuations can be 
read off Fig. [21 However the magnitude of the signal, and 
hence the amplitude of the 2k f charge modulation is sub- 
stantially smaller than in the spin sector. The Luttinger 
liquid parameter, K n of the Hubbard model is bounded 
by 1/2 < K p < 1 [l£|. Hence 4k f charge-charge corre- 
lations which decay as r~ 4Kp are very much suppressed 
in comparison to 2k f charge fluctuations which decay as 
r ~ 1 ~ K p . This stands in accordance with the data of Fig. 
and no divergence in the 4k f charge susceptibility is 
expected. 

Having pinned down energy scales we now consider the 
dynamical spin and charge structure factors: 



S % (k,u) 



1 



-(3uj 



ImXc(fc,o;) 



(4) 



Fig. IHlplots the dynamical spin structure factor as a func- 
tion of temperature. As apparent below the crossover 
scale Tj the two spinon continuum of excitations with 
gapless excitations at 2kp is clearly visible. Furthermore, 
below Tj, a well defined spin velocity can be read off the 
data yielding v s /t ~ 1. This result compares favorably 
with the zero temperature results of Ref. [ll[. Hence, 
in the spin sector Tj marks the temperature scale below 
which the the overall features of the zero temperature 
dynamical spin structure factor become apparent. 

The dynamical charge structure factor is plotted as a 
function of temperature in Fig. 21 Again, below Tj, one 
can read off the charge velocity, v c ~ 1.9 which favorably 
compares with the zero temperature data of [ll[. Given 
the small amplitude of the 2k f charge fluctuations, we 
are unable to reliably pin down the expected gapless ex- 
citations at 2k f as well as at 4k f. 



B. Single particle excitation spectrum 

Our major interest here is to study the temperature 
dependence of the single particle spectral function and 




FIG. 5: Single particle excitation spectrum for T = shown 
as gray scale plot with a logarithmic intensity scale. These 
data stem from the DDMRG calculations of H. Benthien et 
al. 



compare to the experiments of Ref. Q. At the lowest 
temperatures considered in Ref. 0, the photoemission 
results compare favorably with T = DDMRG calcula- 
tions of Ref. 3]. In the photoemission spectra one can 
identify a spinon branch, a holon branch as well a holon 
shadow band. Those features compare well with the zero 
temperature DDMRG data shown in Fig. [5J Let us con- 
centrate on uj/t < relevant for comparison with pho- 
toemission. In the vicinity of the Fermi wave vector and 
at low energies one clearly observes two features (branch- 
cuts) dispersing linearly with velocities v s (spinon) and 
v c (holon). Those velocities stand in good agreement 
with those determined by our analysis of the spin and 
charge dynamical structure factors. Furthermore, and 
at low energies one can identify a feature at 3k f which 
merges at k = with the holon branch. Following Ref. 
|3 one notes that this feature has the same dispersion 
relation as the holon branch but shifted by 2k f. Hence 
the interpretation of a shadow holon branch which stems 
from a holon scattering off a 2k f spin excitation. This 
is very reminiscent of shadow bands in spin-density wave 
approaches of antiferromagnetic Mott insulators. 

The finite temperature spectra we have obtained with 
the QMC are presented in Fig. The question we wish 
to address is at which temperature scale do the features 
of the T = data become apparent. One can observe 
a clear spinon branch for the different temperatures ex- 
cept at (3t = 2 where spectrum is washed out. The holon 
shadow band can also be identified at least in the region 
of wave vectors where < k < kp. As can be seen in 
the DDMRG T = spectrum, the intensity of the holon 
shadow band rapidly decreases for larger wave vectors 
making it very hard to retrace this feature in our spec- 
tra. Another difficulty arises regarding the holon branch 
in our QMC spectra. To analyze our data we use the 
stochastic MEM, with its well-known difficulties to re- 




FIG. 6: Single particle excitation spectrum as a function of 
temperature on a logarithmic intensity plot. From top to 
bottom: fit = 15, 10, 7, 2. 
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FIG. 7: Temperature dependence of spectral functions at (a) 
k = and (c) k = kf. Integrated spectral weight for (b) k = 
and (c) k — kf. 



solve two peaks close in energy. Nevertheless, at our 
lowest presented temperature (3t = 15, a holon branch 
can be identified. Hence, the finite temperature results 
stand in agreement with the statement that below the 
spin scale Tj, the gross features of the zero temperature 
results are apparent. 

The ARPES measurements of Ref. [5| point towards 
spectral weight transfer in the temperature range 60K < 
T < 260K. At the Fermi wave vector and as a func- 
tion of decreasing temperature spectral weight is trans- 
fered from higher (~ 0.7eV) to lower (~ O.leV) excita- 
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tion energies. The numerical simulations do indeed show 
spectral weight transfer however at a temperature scale 
T > Tj. Figs. plots A{k,uj) as a function of temper- 
ature for k = and k = kf. Upon inspection of the 
data, one observes that at high temperatures (/3t — 2), 
the spectral weight is dominantly located at frequency of 
the holon, u/t ~ — 1.5 for k = and u ~0 a,t k = kf. As 
temperature is lowered thereby generating short ranged 
2k f spin fluctuations, spectral weight is shifted over an 
energy scale set by t to form the spin related features. 
At k = 0, this corresponds to the spinon at u/t ~ —0.5 
and at k = kf to the holon-shadow band at u/t c± —2. 
For T < Tj and within the limitations of the stochastic 
analytical continuation, the data shows no further shift 
of spectral weight. 

IV. CONCLUSION 

We have computed the temperature dependence of the 
single particle spectral function for the one-dimensional 
Hubbard model, for a parameter range which has been 
proposed for the modeling of the TCNQ band in TTF- 
TCNQ organics; U/t = 4.9, t = OAeV and n = 0.59. 
This parameter set reproduces well the overall features of 
the photoemission spectra at T = 60K above the Peierls 
temperature. For this parameter set we have identified 
a magnetic energy scale Tj ~ O.lt below which 2k f spin 
fluctuations are enhanced as a function of decreasing tem- 
perature. For T < Tj the overall features of the T = 
spectral function are apparent. In particular no shift in 
spectral weight between holon, holon-shadow and spinon 
branches is observed below this temperature scale. On 
the other hand, for T > Tj spectral weight transfer over 
an energy scale set by t is observed. Very similar con- 
clusions have been reached for the half-filled Hubbard 
model [l2|. With t = OAeV we obtain a magnetic scale 
Tj ~ 400K and hence we are unable to account for the 
spectral weight transfer observed in the photoemission 
experiments in the temperature range 60K < T < 260K 

i. 

Assuming that a pure electronic model is valid to ac- 
count for the temperature dependence of the spectral 



function, other parameter sets are required to under- 
stand the experimental data. The aim is to keep the low 
temperature spectral function similar to that observed in 
this work since it compares well with the low tempera- 
ture photoemission data, but to reduce the spin scale Tj. 
X-ray scattering experiments of Ref. 13] suggest that 
both 4k f and to 2k f charge fluctuations are present at 
low temperatures and that above 150 K only 4k f scatter- 
ing is present. To model such dominant 4k f fluctuations 
one requires a Luttinger liquid parameter K p < 1/2 [H| • 
Since the Hubbard model has 1/2 < K p < 1 additional 
terms such as a nearest neighbor Coulomb repulsion V 
is required. However, we expect that U-terms in the 
Hamiltonian will enhance the overall low temperature 
band width. This band- width problem could be corrected 
by reducing the value of the hopping matrix element to 
it's bulk value, t ~ 0.2eU, as inferred from DFT calcu- 
lations p . In turn this would enhance the value of U/t 
and hence reduce the value of Tj. Further simulations 
are required to confirm this point of view. 

The issue of coupling to the lattice is still open. In par- 
ticular since the system is close to a Peierls transition, 
it is not clear that phonons can be omitted. Further- 
more the line-shapes of model calculations at low tem- 
perature are much sharper than the experimentally ob- 
served. Coupling to the lattice, could account for this 
broadening. 
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